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ABSTRACT 

We present the results of the first survey and monitoring study of the linear polari- 
sation properties of compact, flat-spectrum radio sources at millimetre/submillimctrc 
wavelengths and discuss the implications of the inferred magnetic held structure for 
the emission models involving shock waves in relativistic jets on sub-parsec scales. 
We find significant polarisation in most sources but, in general, the magnetic field on 
sub-parsec scales is less well ordered than on parsec scales. We observe no difference in 
polarisation properties between the BL Lac objects and compact flat-spectrum quasars 
at these wavelengths. Although we find the behaviour of some sources, particularly the 
most highly polarised, to be very consistent with the predictions of transverse shock- 
in-jet models, the detailed behaviour of most objects is not. Conical shock structures 
can more readily explain the observed diverse behaviour of the sample, although some 
degree of bending of the jet may still be necessary in some cases. 

Key words: BL Lacertae objects : polarisation — quasars : polarisation — galaxies 
: jets : shocks — radio continuum : galaxies — submillimetre : JCMT : UKT14 
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Compact, flat-spectrum, radio-loud extragalactic sources, 
sometimes also known as 'blazars' (Angel & Stockman 1980, 
Antonucci 1993), exhibit many common characteristics. In 
particular, they show core-jet morphology in VLBI im- 
ages with a succession of synchrotron self-absorbed knot- 
like structures combining to give the overall flatness of 
the radio spectrum, and frequently exhibit apparent su- 
perluminal motion along with large-amplitude variability at 
all wavelengths from radio to gamma-rays. Significant and 
variable linear polarisation is also seen at all wavelengths 
where measurements have been attempted. At millime- 
tre/submillimetre wavelengths the emission arises through 
the synchrotron mechanism in the core regions of the jets, 
i.e. close to their origin (see e.g. Gear 1988), which are 
completely unresolved on even the highest-frequency VLBI 
maps. Photometric monitoring in this band has greatly im- 
proved our understanding of the physical processes in the 
sub-parsec scale jets and has led to the development of the 
canonical shock-in-jet paradigm (Marscher & Gear 1985 and 
Hughes, Aller & Alter 1985, 1989a, 1989b). An essential in- 
gredient of jet models that has so far been missing is the 



geometry of the magnetic field on these very small scales 
and how it may vary on different physical scales. The mil- 
limetre/submillimetre emission is expected to be linearly po- 
larised if the magnetic field is ordered and, with very rare ex- 
ception (e.g. Robson et al. 1983), it is always optically thin. 
This means that internal synchrotron self-asorption and ex- 
ternal Faraday rotation of the polarisation vector can be 
ignored. Therefore, we can reliably assume that the mag- 
netic field direction must lie perpendicular to the observed 
polarisation position angle (see e.g. Pacholczyk 1970). 

The advent of polarisation-sensitive VLBI imaging has 
been a powerful tool in investigating magnetic field be- 
haviour (e.g. Cawthorne et al. 1993, Gabuzda et al. 1994). 
However, no information has been available on the degree of 
order and the orientation of the magnetic field close to the 
point of creation of the jet. In the early 1990s we began a 
1.1 and 0.8 mm polarisation survey of a significant sample 
of flat-spectrum radio sources, most of which have subse- 
quently been measured over several epochs at 1.1 mm. In 
this paper we present the data from this survey, which allow 
us to carry out the first investigation of the magnetic field 
properties in the core of compact radio jets, and we discuss 
the implications for models of these sources. 
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2 SAMPLE 

Polarisation detections have been obtained on 26 different 
sources, half of which have been observed at 4 or more differ- 
ent epochs between 1991 and 1996 and single-epoch upper 
limits have been measured on a further two objects. A com- 
plete list of the objects with their polarimetric observations 
to date is given in Table 1. The main criterion adopted for 
the selection of sources has been that their flux density at 
1.1 mm should be at least 1 Jy so that a detection can be 
achieved in a reasonable time. Due to the highly variable na- 
ture of these objects and in order to ensure that changes in 
the emission between epochs are detected above the uncer- 
tainty level of the measurements, a 5a polarimetric detection 
was always attempted and often achieved. Other selection 
criteria include the availability of published high-frequency 
VLB1 observations, access to the target sources at various 
times of the year and keeping a balanced mixture of different 
optical types. 

3 OBSERVATIONS 

3.1 Data Acquisition 

The observations were acquired using the Aberdeen/QMW 
Polarimeter (Murray 1991 and Murray et al. 1992) in con- 
junction with the UKT14 bolometric receiver (Duncan et al. 
1990) located at the east Nasmyth focus of the 15-m James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea, Hawaii. 
The polarimeter consists of a rotating half-waveplate fol- 
lowed by a fixed wire-grid analyser, ft is operated in the 
'Step and Integrate' mode, whereby a standard photometric 
measurement is performed at a number of equi-spaced posi- 
tions of the waveplate (between 0° and 360°), thus modulat- 
ing the polarised component in the signal. The observations 
are automated by adopting a waveplate rotation and beam 
switching scheme which drives the waveplate and nods the 
telescope at each waveplate position. The telescope's sec- 
ondary mirror is used to chop between the source and the 
sky thermal background at the default frequency of 7.8125 
Hz with a chop throw of 60". The position angle of the lin- 
ear polarisation vector is rotated by the parallactic angle of 
the source due to the alt-azimuth mounting of the telescope 
and it is further rotated by its elevation in being reflected 
off the tertiary mirror. These rotations can be corrected for 
either during acquisition or off-line in software. When mak- 
ing purely photometric observations with the polarimeter 
placed in front of the receiver, two separate measurements 
with the waveplate at 0° and 45° are taken and averaged. 
Further technical details are given in Murray (1991) and 
Nartallo (1995). 

3.2 Data Reduction 

The raw output of the receiver S can be represented as a 
function of the waveplate position 9 by the following equa- 
tion (e.g. Nartallo 1995), 

S(0) = i[/ u +/ P [l+cos(40-2 X )]] (1) 

where I u and I p are the intensities of the unpolarised and 
polarised components of the signal respectively, and x is 
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Waveplate Position 8 (°) 

Figure 1. Examples of the fitting and noise removal procedures. 
Plot (a) shows a good fit to a noise-less first half-cycle of the 
waveplate, while plots (b) and (c) show fits to the second half- 
cycle of the same observation before and after applying the noise 
removal algorithm. 

the position angle of the linear polarisation vector. In terms 
of the Stokes parameters I, Q and U (which fully describe 
linear polarisation) this can be written as: 

S(0) = i[J + Qcos(40) + Usin(40)] (2) 

where I = I u + I p , Q = J p cos(2x) and U = 7 p sin(2x)- 
The Stokes parameters are estimated by performing a least 
squares fit on a number of data points taken at different po- 
sitions of the waveplate. This technique allows the removal 
of individual spikes which are occasionally observed in the 
raw data due to changing sky emissivity, cosmic rays or elec- 
tronics glitches. A noise removal algorithm was developed 
that optimises the \ 2 value of an n-point fit while maxi- 
mizing the number of points used within a waveplate cycle. 
Other methods tried for determining the Stokes parameters 
include fitting the data by a Fast Fourier Transform and 
the estimation of Q and U directly from photometric mea- 
surements taken 45° apart. However, spiked points cannot 
be taken out when using these techniques, often resulting in 
unusable estimates of Q and U. 
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In practice, a typical observation consists of 16 wave- 
plate positions that are fitted in half-cycles of 8 points. Half- 
cycles that have more than 3 points removed are thrown 
away. Figure 1 illustrates a fit to a 16-point observation of 
a highly polarised source (~ 15% with position angle 40°). 
Plot (a) shows the normalised raw data for the first half- 
cycle and its corresponding fit and plot (b) is the equivalent 
representation for the second half-cycle, which is clearly af- 
fected by noise. The peak of the fitted curve gives the per- 
centage polarisation and its displacement from the 8 — 0° 
position is half the position angle of the polarisation vector. 
While the fit shown in plot (a) gives a degree of polarisa- 
tion of 13% and a position angle of 45°, the values obtained 
from the second half-cycle shown in plot (b) are 31% and 
61° respectively. Plot (c) shows the result of applying the 
noise removal algorithm to (b) and it is now seen that it 
is in good agreement with (a) (giving 15% polarisation and 
38°). 

After successfully fitting the raw data the normalised 
Stokes parameters q and u are calculated and corrected for 
the Instrumental Polarisation (IP) and rotations. The main 
contributor to the IP is the protective wind blind of the 
JCMT (Flett & Murray 1991) and it is therefore dependent 
on the actual blind used and its tension at each epoch (fur- 
ther details of the origin and characteristics of the IP are 
discussed in Murray 1991). The value of the IP to be sub- 
tracted from observations of point sources such as blazars 
is estimated from observations of planets with small angu- 
lar size (e.g. Uranus, Mars or Neptune) which are consid- 
ered to be unpolarised at millimetre wavelengths. Typically, 
the IP measured on these planets has values in the range 
1.5 — 1.8% at 1.1 mm and close to 0.5% at 0.8 mm, which 
remain largely unchanged over the period of an observing 
run. Unlike the IP position angle (which tracks with eleva- 
tion) the level of linear instrumental polarisation does not 
have an elevation dependance. The IP position angle at zero 
elevation is ~ 173°. 

Independent estimates of q and u are obtained from sev- 
eral revolutions of the waveplate and checked for homogene- 
ity by performing a Kolmogorov-Smirnov test, that com- 
pares their distributions with a gaussian distribution and 
can indicate the presence of spurious observations. All the 
observations that are finally accepted are averaged to im- 
prove the statistical determination of the polarisation pa- 
rameters. The number of observations available for each 
source and their consistency allow an error weighted average 
to be used in most cases, which improves the final signal to 
noise. If only a few measurements are available and/or their 
individual signal to noise varies significantly between them, 
a straight mean is used instead. 

Finally, the degree of linear polarisation P and its un- 
certainty a p can be obtained from the averaged normalised 
Stokes parameters using the expressions, 



Pobs = \/q 2 + m 2 

y/(q X a q ) 2 + (u X a u ) 2 



obs 



(3) 



(4) 



the observed degree of polarisation P b s is then corrected 
for bias using the Wardle and Kronberg method (Wardle & 
Kronberg 1974) which is adequate for P b s /u p > 0.7 (Sim- 
mons & Stewart 1985), 



P = y/PL-°l (5) 

The position angle \ and its uncertainty a x are given by, 
X= \ arctan (^j (6) 

a *=2p- (?) 
where the expression for a x is a good approximation if P 3> 
cr p (Naghizadeh-Khouei & Clarke 1993). (A more detailed 
account of all the data reduction algorithms used and tests 
carried out on them is given in Nartallo 1995). 



3.3 Observing Wavelengths 

Polarimetric observations at 1.1 and 0.8 mm have been ob- 
tained at nine and three different epochs respectively. A 
few extra observations of selected sources taken during the 
two commissioning runs of the polarimeter (Flett & Mur- 
ray 1991) have been added to the database and re-analysed. 
Initially, observations were systematically made at the two 
wavelengths, thus limitting the number of objects that could 
be covered in a single epoch. Data from these earlier epochs 
show that the polarisation properties of most sources at 1.1 
and 0.8 mm are very similar. This is illustrated by Figure 2, 
where the variability with epoch of the percentage linear 
polarisation and position angle is shown for three sources 
at these two wavelengths. The plots include four-epoch ob- 
servations of 3C273 and three-epoch observations of 3C279 
and 3C454.3. These three objects have different overall po- 
larisation properties, ranging from lowly polarised (3C273) 
to very highly polarised (3C279) . Since there is little knowl- 
edge to gain by monitoring the polarisation at the two wave- 
lengths, it was decided to concentrate on single-frequency 
polarimetry at 1.1 mm in order to make the sample statisti- 
cally significant by increasing the number of sources measur- 
able during an observing run. The choice of wavelength was 
determined by the fact that a good signal to noise can be 
achieved more quickly at 1.1 mm and also this atmospheric 
window is less affected by weather conditions. However, the 
level of instrumental polarisation is a factor 3 greater than 
at 0.8 mm. 

Whenever possible, photometric measurements of the 
sources have also been taken at 2.0, 1.1 and 0.8 mm, from 
which a three-point spectral index is calculated (the source 
spectrum I(y) oc v a ). For the earlier epochs, the photom- 
etry has been extracted from the 1.1 and 0.8 mm polari- 
sation observations and the spectral index is therefore es- 
timated from two wavelengths only (thus having greater 
uncertainties). The photometric measurements are flux cali- 
brated against planets measured close in time and at similar 
airmasses (Griffin el al. 1986, Orton et al. 1986, Griffin & 
Orton 1993). 



4 RESULTS 

The polarisation and photometric measurements obtained 
on each source during the course of this monitoring pro- 
gramme are tabulated in Table 1 for each observation epoch. 
The table includes the percentage linear polarisation P and 
position angle \ at 1.1 and 0.8 mm, the flux density F in 
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Figure 2. Comparison between the 1.1 and 0.8 mm polarisation properties of three compact, fiat-spectrum radio sources. The data 
represented by triangles joined by a solid line correspond to observations of 3C273, the stars and dashed lines correspond to 3C279 and 
the squares and dashed-dotted lines arc measurements of 3C454.3. The error bars shown represent the uncertainty in each individual 
detection, as given in Table 1. 



Janskys at 2.0, 1.1 and 0.8 mm and the spectral index a 
determined from the available wavelengths. In Figure 3 we 
have plotted the four measured parameters against epoch 
of observation for every source detected. The polarisation 
position angle plots are centred about the position angle of 
the VLBI jet as given in Table 2, indicated by a horizontal 
dashed-dotted line; the direction perpendicular to the jet is 
indicated by two dotted lines at 90° . The inferred magnetic 
field geometry is perpendicular to the jet if x lies close to 
the VLBI jet direction and parallel if they are closer to the 
dotted lines. A summary of the main features observed in 
the 1.1 mm measurements obtained on each source is given 
in Appendix A. 



4.1 Measurements of Spectral Index as a Guide 
to Opacity 

Photometric data acquired at three millimetre wavelengths 
have been used to estimate the spectral index. A flat spectral 
index suggests that the emission is optically thick, especially 
if it coincides with a significant increase in the total flux. If 



there is significant opacity, synchrotron self-absorption can 
swing the polarisation position angle by 90° (Pacholczyk 
1970). There are 12 sources with at least three epochs of 
simultaneous measurements of a, F, P and x (see data in 
Table 1 and Figure 3) which allow us to test for opacity. 

The variability in the spectral index measurements 
is different between sources. In objects 0235+164, OJ287, 
1749+096 and 1921—293, the spectral index varies within a 
small range of values —0.7 < a < —0.3, despite significant 
changes in F, P and x- in 3C454.3, however, the spectral 
index takes values — 1.8 < a < —0.2, and although position 
angle swings by as much as 90° have also been observed, nei- 
ther of these changes can be associated with 'flares' in the 
source emission, as the flux remained constant over the first 
five epochs and the polarisation level varied only slightly. 
We do find a tendency for a flatter spectral index as the flux 
increases in 3C279, 1335—127 and BL Lacertae, but while 
the polarisation also increases in 3C279, it decreases in the 
other two sources. In all three objects, however, x remains 
largely unchanged, which is inconsistent with the emission 
becoming optically thick during a flare. In 0420—014 the 
epoch with flatest spectral index coincides with the highest 
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Table 1 . Polarimetric and photometric data tabulated in terms of bias-corrected percentage polarisation P, position angle \ and flux density F. The subscript indicates the wavelength 
in mm. The epoch is expressed in decimal years. The spectral index a of the mm source spectrum is calculated as a 'positive' power of the frequency oc v OL ). 
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.4 
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1993.15 
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2.4 


8 


7 






1.00 


0.10 










1335-127 


1992.81 


5.6 


0.6 


116 


3 






2.50 


0.23 












1993.15 


6.2 


0.7 


120 


3 


4.90 


0.57 


3.10 


0.32 


2.56 


0.27 


-0.71 


0.17 




1994.15 


2.5 


0.7 


122 


8 


6.02 


0.57 


4.12 


0.11 


5.08 


0.32 


-0.14 


0.12 




1994.40 


7.9 


2.3 


121 


8 


4.98 


0.38 


3.23 


0.25 


2.24 


0.19 


-0.85 


0.12 




1995.90 


6.0 


1.1 


139 


5 






4.93 


0.50 
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1992.02 


3.7 


1.1 


33 
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1.59 


0.16 












1995.90 


9.1 


3.9 


13 


11 






1.10 


0.22 
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7.8 


2.2 
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0.14 
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5.5 


1.6 


4 
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1633+382 
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7.9 
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0.2 
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2 
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47 
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2.5 


0.4 
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4 


2.8 0.7 
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5.2 


0.5 


57 
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1.0 


95 
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1.3 


59 


7 








1994.40 


11.2 


2.9 


61 
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1995.98 


4.5 


1.2 
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8 






1741-038 


1993.15 


4.6 


1.6 
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9 
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3.7 


0.4 


10 
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1993.15 


4.7 


0.7 
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8.5 


2.8 


62 


9 
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1992.81 


12.8 


1.3 


40 
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1991.43 
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0.3 
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1 


10.2 0.6 


103 




1993.15 


9.7 


0.6 
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1994.15 


6.9 


1.0 
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1994.40 


10.4 


1.8 


134 


5 
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0.7 


23 
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8.1 
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168 
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0.9 
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7.2 


1.3 


20 
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1994.15 


11.8 


3.6 


19 


8 








1995.59 


3.0 


0.7 


17 


7 
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1990.42 


3.8 


0.9 


27 


6 








1992.81 


2.2 


0.9 
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11 
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1990.42 


4.5 


0.3 


109 


2 


6.1 0.9 


110 




1991.43 


2.9 


0.3 


15 


3 


2.0 1.3 


28 




1992.16 
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1.6 


47 


6 


10.2 0.6 
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0.4 
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0.16 
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Figure 3. The mm spectral index and the flux, percentage polarisation and position angle measured at 1.1 mm is plotted against epoch. 
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Figure 3 - continued 
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Figure 3 - continued 
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Figure 3 - continued 
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Figure 3 - continued 



flux and polarisation measured on this source. In 0735+178, 
on the other hand, the lowest a is measured when F and P 
are both low. In other sources the behaviour is even more 
complex. Take for example the last three epochs of 3C345 
with spectral index information, obtained in just over a year 
(between February 1993 and May 1994): the flux, that had 
been dropping steadily since 1991, falls from 3.2 to 1.7 Jy, 
yet a changes from —0.82 to —0.14 to —1.15 and P increases 
from 2.9 to 5.4 to 11.2%. These drastic changes in spectral 
index do not seem to be directly related to changes in flux. 
3C273 also exhibits complex behaviour and, although half of 
the spectral index values available were obtained from two 
wavelength measurements and therefore have large uncer- 
tainties, we have good signal to noise in all F and P mea- 
surements. We have data covering two flaring events during 
which the polarisation was high. The epoch of lowest polar- 
isation (0.7%) and largest position angle swing (60°), coin- 
cides with the lowest flux measured (11 Jy) and a steep spec- 
tral index of —0.94 ±0.45. The measurements of the less fre- 
quently observed sources in the sample all give equally steep 
spectral indices, with the possible exception of 1514—241 
(q = 0.05). 

None of this evidence is suggestive of optically thick 
emission and we therefore make the asumption that all 
sources have optically thin synchrotron spectra. 



4.2 Degree of Polarisation 

Examination of the database shows that the degree of po- 
larisation P can vary considerably in most sources, even in 
timescales of a few weeks (see e.g. the change in the 1.1 
mm percentage polarisation in sources 0420—014, 0735+178, 
OJ287, 3C273 & 3C345 between epochs 1992.02 and 1992.16 
given in Table 1). However, there seems to be a characteris- 
tic range of polarisation values for each source: for example 
in 3C273 0.7% < P < 6.8%, in 3C279 7% < P < 17% and in 
3C454.3 3% < P < 8.5. We have estimated the mean polar- 
isation P for 20 different sources by averaging the measure- 
ments over all epochs available for each object. There are 
another 6 sources measured only once, and 2 further objects 
on which we have la upper limits. The polarisation values 
obtained are listed in Table 2. 

In particular, we find that the millimetre emission from 
OJ287, 3C279, 1921-293 and BL Lacertae is always highly 
polarised, with P in the range 9 — 11% in all cases, while 
that from 3C273 and 3C84 typically remains at much lower 
levels (2 — 3% or less). For a majority of sources, though, the 
mean polarisation level is typically 4 — 7%. Figure 4 shows 
the distribution of P for multi-epoch sources (dark bars), to 
which we have added the measurements from single-epoch 
sources (shaded bars). The main peak of this distribution is 
~ 4.5% and there is a strong secondary peak at ~ 10%. 
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Figure 4. Distribution of percentage polarisation plotted in 3% 
bins. The dark bars represent the mean polarisation values ob- 
tained on sources with multiple observations and the light bars 
show the data on single-epoch sources. 



We show below that these two broad groups of sources 
with 'high' and 'intermediate' polarisation levels behave 
somewhat differently in terms of the changes in the polari- 
sation position angles measured. These data also show that 
although the magnetic fields are certainly ordered on sub- 
parsec scales, where the millimetre emission is coming from, 
the degree of ordering is much less than further out in the 
jets, where polarisation-sensitive VLBI has revealed regions 
with polarisation levels close to the theoretical maximum 
for synchrotron radiation (~ 70%). It is also possible, how- 
ever, that the fields are ordered on even smaller scales and 
that the polarisation detected is the average of several differ- 
ent regions. Only polarisation-sensitive VLBI at 1mm would 
unambiguously settle this question. 



4.3 Relationship between Degree of Polarisation 
and Total Flux 

It is possible to investigate the relationship between the de- 
gree of polarisation and total flux of those sources that have 
been observed at a number of different epochs. Figure 5 is 
a plot of these two parameters for all sources with at least 
three measurements. For some of the most frequently ob- 
served, namely 3C273, 3C279 and 3C454.3, a strong and 
statistically significant positive linear correlation is found 
between the polarised and total intensities (with r > 0.8 and 
confidence level > 99.4 per cent). Very tentative positive cor- 
relations between P and F may also occur in 0420—014 and 
0735 + 178. In other sources, however, there is some evidence 
that these parameters may be anti-correlated, for example 
in 0235+164, 3C345, 1749 + 096 and BL Lacertae. OJ287 is 
a special case since, although no statistical correlation can 
be found when all 8 observations available are considered, 
this is due to a change in the sense of the correlation from 
positive to negative at an epoch when the flux density is 
a maximum and the polarisation is at minimum level (see 
Appendix A for further details). In 4C39.25, 1334-127 and 
1921—293, however, there is no clear relationship between 
these two parameters. 
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Figure 5. Percentage polarisation is plotted against flux density 
for all sources on which there are at least three epochs of data. 
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Table 2. Source name, optical type, mean percentage polarisation and VLBI information is given for all sources detected. Column 1 
lists the sources observed, some with alternative names commonly used. Columns 2 and 3 give the optical identification and references 
in terms of the following classes: BL Lac, HPQ (High Polarisation Quasar), LPQ (Low Polarisation Quasar) and RG (Radio Galaxy). 
Columns 4 and 5 are the mean percentage polarisation at 1.1 mm and standard deviation for sources with multi-epoch observations. 
There are eight single-epoch sources (two of which are la upper limits) for which the polarisation measured and its error are quoted. 
The number of polarimetric epochs available on each source is given in column 6. Columns 7—11 give the position angle of the jet as 
determined from VLBI maps, the distance from the core, frequency of the VLBI observation, references and comments. 



Source Name Optical ID 







Class 


Ref. 


0133+476 




HPQ 


[cl 

L L J 


0235+164 




BLLac 


fbl [el 

L^l i L&J 


0316+413 


(3C84) 


RG 


fbl, [el 

i L&J 


0415+379 


(3C111) 


RG 




0420-014 




HPQ 


fcl fgl 

L^J i L&J 


0430+052 


(3C120) 


RG 


[al [el 

l a \ i L&J 


0528+134 




LPQ 


[cl 

L^J 


0735+178 




BLLac 


fbl, fgl 

L^J i L&J 


0736+017 




HPQ 


fcl, [el 

L^J i L^J 


0829+046 




BLLac 


faljbl 


0851+202 


(OJ287) 


BLLac 


[b],[g] 


0923+392 


(4C39.25) 


LPQ 


[c] 


1055+018 




HPQ 


[b],[c] 


1226+023 


(3C273) 


LPQ 


[b],[c] 


1253-055 


(3C279) 


HPQ 


[e],[c] 


1308+326 




BLLac 


[b],[d] 


1335-127 




HPQ 


[b],[c] 


1413+134 




BLLac 


[a] 


1514-241 




BLLac 


[a] 


1633+382 




LPQ 


[c],[g] 


1641+399 


(3C345) 


HPQ 


[c],[e] 


1741-038 




HPQ 


[c] 


1749+096 




BLLac 


[a],[b] 


1823+568 




BLLac 


[d] 


1921-293 




HPQ 


[b],[g] 


2200+420 


(BL Lac) 


BLLac 


[b],[g] 


2223-052 


(3C446) 


BLLac 


[f],[g] 


2251+158 


(3C454.3) 


HPQ 


[c],[c] 



Polarisation 


VLBI 


T mas 


GHz 


% 


(J 


n 








10.5 


1.8 


1 


165 


1.3 


5 


4.9 


1.0 


4 


-150 


0.5 


22 


1.4 


0.5 


2 


-150 


< 1.0 


43, 100 


<0.7 




1 








5.0 


0.6 


4 


-150 


< 2.0 


43 


<2.2 




1 


-100 


< 20.0 


5 


3.9 


1.0 


2 


70 


0.3 


22 


6.3 


1.3 


4 


73 


0.4 


22 


9.8 


4.6 


2 


-80 




22 


3.9 


1.0 


2 








10.5 


1.4 


8 


-120 


< 0.1 


43, 100 


4.7 


1.4 


4 


—88 


2.2 


43 


9.1 


1.0 


2 


-55 


8.0 


1.7 


2.1 


0.4 


10 


-110 


< 0.3 


43, 100 


10.5 


1.0 


11 


-135 


< 0.1 


100 


9.9 


2.4 


1 


-46 


< 0.7 


5 


5.5 


1.3 


5 


135 


< 0.2 


5 


6.9 


2.3 


3 


-100 


< 8.0 


8.4 


5.5 


1.6 


1 








7.9 


2.7 


1 


-45 


1.3 


5 


4.8 


1.4 


8 


-55 


< 0.2 


43, 100 


4.6 


1.6 


1 


-155 




22 


5.6 


1.8 


3 


65 


~ 1.0 


5 


12.8 


1.3 


1 


-166 


0.8 


5 


9.1 


0.9 


4 


-40 




43 


9.1 


1.2 


6 


-160 


< 1.5 


100 


3.0 


1.0 


2 


-142 


~ 0.1 


100 


4.6 


0.8 


6 


-65 


0.6 


43 



Refs. 


Comments 


fl9l 


Low resolution map 


161 


Weak component 


[21,191 


High resolution maps 




No VLBI maps available 
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Same PA at lower frequencies 
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No other maps available 




No VLBI maps available 
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PA similar to that at 5 GHz 


[i] 


Unresolved core 


[13] 


Low frequency maps 
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Several maps available 


[2] 


Well determined PA 


[5] 


Single VLBI observation 


[18] 


Core unresolved 


[14] 


Low resolution 




No VLBI maps available 


[4] 


PA of inner jet 


[2] [10] [16] 


Large range of PAs 


[12] 


Short, weak jet 


[18] 


Unresolved, East elongation 


[4] 


Agrees with earlier maps 


[17] 


No other maps available 
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Same PA at lower frequencies 


[11] 


PA of innermost jet 


[7] 


Same PA at lower frequencies 



OPTICAL ID REFERENCES - [a] Burbidge & Hewitt (1992); [b] Ghisellini et al. (1993); [c] Impey & Tapia (1990); [d] Kiihr & Schmidt (1990); [e] Moore 
& Stockman (1984); [f] Rusk (1990); [g] Wchrlc et al. (1992). 

VLBI REFERENCES - [1] Albcrdi et al. (1993); [2] Baath et al. (1992); [3] Bloom et al. (1996); [4] Cawthornc et al. (1993); [5] Gabuzda et al. (1992); [6] 
Jones et al. (1984); [7] Kcmball et al. (1996); [8] Krichbaum et al. (1990); [9] Krichbaum et al. (1992); [10] Krichbaum et al. (1993); [11] Lcrncr et al. (1993); 
[12] Marschcr et al. (1997); [13] Padriclli et al. (1986); [14] Pcrlman et al. (1994); [15] Pilbratt et al. (1987); [16] Rantakyro et al. (1992); [17] Shcn et al. 
(1997); [18] Wchrlc et al. (1992); [19] Zcnsus et al. (1984); [20] Zhang & Baath (1991); [21] Zhang et al. (1994). 



These results can shed light into the physical processes 
taking place in the most compact regions of the jets and 
could also provide clues about their morphology, as will be 
discussed later in §5. 

4.4 Magnetic Field Orientation 

In order to relate these polarisation observations to the mag- 
netic field geometry in the source, we make the assumption 
that the emission is optically thin synchrotron radiation, 
for which the polarisation position angle is perpendicular 
to the magnetic field direction in the source frame. The 
magnetic field geometry can then be related to the posi- 
tion angle of the jet axis, as determined from VLBI observa- 
tions. Many core dominated sources, however, show signifi- 
cant bending from parsec to kiloparsec scales (e.g. see 3C84, 



3C273, 3C279, 3C345, 3C454.3 & 4C39.25 in Readhead et 
al. 1983) and we know that if further bending occurs be- 
tween the scale of the centimetre emission and the millime- 
tre emission, our conclusions could be severely compromised 
or even completely erroneous. Even for sources with appar- 
ently straight jets from parsec to kiloparsec scales, we cannot 
rule out the possibility of bending on smaller scales. Ideally 
we would like to compare our polarisation data with VLBI 
maps obtained at 1 mm, but in practice we are foced to 
estimate the direction of the jet structure from the highest- 
frequency VLBI map available to us. These have been listed 
in Table 2, where we have also included the distance from 
the core to the nearest component (in mas), the frequency 
of the observation and a brief comment on the character- 
istics of the map. In the rest of this discussion we refer to 
the magnetic field orientation thus determined, instead of 
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Figure 6. Percentage polarisation as a function of magnetic field orientation with respect to the inner jet direction. Different symbols 
have been used to represent measurements of sources whose mean polarisation P lie in three different ranges: empty stars for P > 9%, 
solid circles for 9% > P > 3% and empty triangles for P < 3%. Measurements of single-epoch sources have also been included. The error 
bars have been omitted for clarity but all observations are better than 3o\ 



the measured polarisation position angle \- The convention 
used throughout is that the magnetic field is perpendicular 
to the VLBI jet when | Qvlbi — X I is 0° an d parallel at 90°. 

In Figure 6 we plot degree of polarisation P against 
| Qvlbi — X • We treat measurements of the same source 
at different epochs as independent points, eventhough some 
sources have been observed more frequently than others. 
This allows a realistic picture of the overall behaviour to 
emerge, rather than taking an average of the positon angle 
measured for a source which might show radically different 
behaviour at different times. We have used, however, differ- 
ent symbols to represent sources whose mean polarisation P 
lie in three different ranges as discussed above: empty stars 
for P > 9%, solid circles for 9% > P > 3% and empty trian- 
gles for P < 3%. Measurements of 8 sources observed only 
once have also been included as belonging to one of these 
three categories. 

This figure clearly shows that the more highly polarised 
sources are rarely found with polarisation levels lower than 
6% and that there is a strong tendency for them to have 
magnetic fields largely perpendicular to the jet axis. Only 
a sixth of the measurements in this subset of objects have 
| Qvlbi — x \> 45°, and 5 of the 'starred' data points close 
to 45° belong to observations of OJ287 alone. 



Sources with an 'intermediate' level of polarisation do 
not have a preferred magnetic field orientation. In fact, many 
of the most frequently observed sources in this category ex- 
hibit magnetic fields either parallel or perperdicular to the 
jet at different epochs, e.g. 0235+164, 0735+178, and any 
orientation in between, such as 4C39.25, 3C345 and 3C454.3 
where the position angle measured is clearly seen to change 
gradually (in tens of degrees) rather than in 90° swings. 
However, the magnetic field in 1335—127 is clearly perpen- 
dicular to the jet and does not change over 5 epochs spanning 
more than 3 years. 

It is also clear from Figure 6 that the measurements 
from sources with low mean polarisation are largely associ- 
ated with magnetic fields closer to the aligned configuration 
rather than perpendicular to the jet. However, this must 
be treated with caution as 10 of these points correspond to 
observations of 3C273 alone. 

All of the examples quoted above can be verified by 
looking at the raw data plotted in Figure 3. One and 
two-dimensional Kolmogorov-Smirnov (K-S) tests of P and 
| Qvlbi — X I between any two categories all give high K-S 
statistics with confidence levels > 99.9 per cent, confirm- 
ing that these groups have statistically distinct polarisation 
properties. 
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4.5 Polarisation Properties of BL Lac Objects 
and Flat-Spectrum Radio Quasars 

4-5.1 Source Classification 

Although all the sources in the sample have so far been gen- 
erally treated as 'blazars', i.e. compact, flat-spectrum radio 
sources with highly variable and polarised emission, most 
objects have traditionally been classified as belonging to one 
of two optical types, namely BL Lac objects (BL Lacs) and 
Flat-Spectrum Radio Quasars (FSRQs). The latter can be 
divided into High Polarisation Quasars (HPQs) (also known 
as Optically Violently Variable quasars or OVVs) and Low 
Polarisation Quasars (LPQs) according to whether their op- 
tical polarisation is > 3% (Moore & Stockman 1981). The 
classical distinction between BL Lacs and HPQs is now un- 
clear as emission lines have been observed in many BL Lacs 
(see e.g. Antonucci 1993 for a discussion). 

In Table 2 we have indicated the most common classifi- 
cation for each source, and we shall adopt this classification 
for the purpose of comparing the properties of these three 
types. Of the 28 sources in our sample we have identified 
11 BL Lacs, 10 HPQs and 3 LPQs. However, some ambi- 
guity exists about the correct classification of some objects, 
especially between the BL Lac and HPQ categories. In fact, 
0133+476 and 1335—127 have at times been classified as BL 
Lacs (Zensus et al. 1984 and Kuhr & Schmidt 1990, respec- 
tively) and 1308+326, 1823+568 and 3C446 have all been 
classified as HPQs in a number of studies (e.g. Moore & 
Stockman 1981, Impey et al. 1991, Burbidge & Hewitt 1992, 
Ghisellini et al. 1993). Also, the BL Lac object 0235 + 164 
has been proposed as a candidate micro-lensed quasar (Saust 
1992). Sources 3C84, 3C111 and 3C120 are usually described 
as radio galaxies. 



4-5.2 Comparison between Different Optical Types 

One of the initial results from the first polarisation-sensitive 
VLBI surveys at centimetre wavelengths was the claim that 
BL Lacs and FSRQs differed in the orientation of their mag- 
netic fields with respect to the jet axis, being always perpen- 
dicular in the BL Lacs and always parallel in the quasars. 
This claim has recently been moderated somewhat, however 
it is still quite marked on parsec scales (see e.g Cawthorne 
et al. 1993 and Gabuzda et al. 1994). Our sample contains 
a similar number of sources of either class (although we 
shall differentiate between the HPQs and the LPQs amongst 
the quasars) and there are multiple observations (obtained 
at different observing epochs) of most objects, making the 
number of measurements in our database statistically sig- 
nificant. We are aware, however, that a few sources have 
been measured more often than others and that this intro- 
duces a certain degree of bias; we have tried to point this 
out wherever it is relevant. 

In Figure 7 we compare the distributions of (a) polari- 
sation level, (b) magnetic field orientation, (c) flux density 
and (d) spectral index for BL Lacs (black bars) , HPQs (dark 
shaded bars) and LPQs (light shaded bars). The bin sizes 
of the histograms are indicated by the horizontal axis of the 
plots: 3% in P, 15° in | 9vlbi — X l> 3 Jy in F and 0.3 in a. 

The shape of the distributions of P and | 9vlbi — X I are 
'visually' very similar between the BL Lacs and the HPQs. 



For both classes, 80 per cent of the polarisation measure- 
ments lie in the range 3% < P < 12% and the magnetic field 
can take any orientation between the perpendicular and the 
parallel to the jet axis. However, while the mesurements of 
BL Lacs are equally distributed between those with mag- 
netic fields oriented between 0° — 45° and those between 
45° — 90°, the HPQs favour the perpendicular configuration 
with a ratio >3:1. In fact, there are no HPQ measurements 
in the 45° — 60° bin, while 5 different sources contribute to 
the 7 measurements of BL Lacs in this bin. We also find con- 
tributions from 6 BL Lacs and 7 HPQs in the 0° — 15° bin, 
and from 6 BL Lacs and 8 HPQs in the 15° - 30° bin. This 
shows that the multiple data points from most sources tend 
to spread throughout the plot. In contrast, there are 9 mea- 
surements of 3C279 (the source most frequently observed) 
in the first bin alone, 5 measurements of OJ287 between 
30° — 60° and 6 measurements of 3C345 between 60° — 90° , 
which influence the concentration of data points in these 
areas. The data available on the LPQs (15 measurements) 
come from 3 different sources, but 3C273 contributes with 
10 entries and 4C39.25 with 4. The distributions of P and 
| Ovlbi — X I f° r this class are therefore strongly biased 
by the characteristics of a single source and thus any com- 
parison with the other classes must be treated with care. 
Although all orientations of the magnetic field have been 
observed at least once in this group, there is a concentra- 
tion of measurements close to direction parallel to the jet. 
Also, half of the polarisation measurements are in the < 3% 
bin and none have been obtained above 8%. 

A K-S test of the P parameter between BL Lacs and 
HPQs gives a K-S statistic of 0.14 with significance level 
of 16 per cent, indicating that the two samples being com- 
pared could have been obtained from a single population of 
objects. A very similar result is obtained on the | Qvlbi — X 
parameter (0.13 with 11 per cent significance) and in a 2- 
dimensional K-S test combining these two parameters (0.16 
with 21 per cent significance). Similar tests carried out be- 
tween either of these two classes and the LPQs all give K-S 
statistics > 0.5 with significance levels of > 99 per cent, 
clearly differentiating between them. 

We therefore find no evidence of any significant differ- 
ence in the properties of the magnetic field (degree of order 
and orientation) between BL Lacs and HPQs on the scales 
sampled by millimetre/submillimetre measurements. How- 
ever, there is some evidence that the magnetic field in the 
inner jets of LPQs may be generally less-ordered and be 
preferentially aligned with the jet, although this needs to be 
confirmed by increasing the number of different LPQs in the 
sample. 

Pannels (c) and (d) of Figure 7 summarise the informa- 
tion obtained from the photometric data. The distribution 
of flux density measurements is very different between BL 
Lacs and HPQs. Practically all the BL Lac measurements 
are <3 Jy, however the flux distribution for the HPQs peaks 
between 3 — 6 Jy and three different sources have been mea- 
sured with F > 6 Jy. The flux distribution of the LPQs is 
difficult to asses as it is dominated by the measurements 
of 3C273. K-S tests between any two categories all give 
a statistic of 0.6 with significance level > 99.9 per cent, 
thus differentiating between the three classes. On the other 
hand, the distribution of spectral index is similar between 
BL Lacs and HPQs and the K-S test fails to separate the 
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Figure 7. Comparison between the distributions of (a) percentage polarisation, (b) magnetic field orientation, (c) flux density and (d) 
spectral index for the BL Lacs (black bars), the HPQs (dark shaded bars) and the LPQs (light shaded bars). Note that the bin sizes in 
all four histograms are the same for all three types being compared and are as follows: (a) 3%, (b) 15°, (c) 3 Jy and (d) 0.3. Also note 
that all the measurements available in the database have been included and therefore most sources have multiple entries. 



three classes (with significance level of < 75 per cent in all 
cases). However, the BL Lacs measurements average to a 
somewhat flatter mean spectral index (—0.56 ± 0.35) than 
the HPQs (-0.70±0.36), which is similar to that of the LPQs 
(—0.72 ± 0.33), and although the difference in the spectral 
slope may not be statistically significant, this result is in 
in very close agreement with that of Gear et al. (f994) ob- 
tained on a larger sample of objects using a greater number 
of wavelengths. 

The relationship between P and F (discussed in §4.3) 
could provide very tentative evidence of different behaviour 
between the BL Lacs and the quasars: in general, P and 
F seem to have a tendency to be anti-correlated in BL 
Lacs (e.g. 0235+164, some epochs of OJ287, 0829+046, 
1413+135, 1749+096 and especially BL Lacertae), whereas 
quasars are more likely to show a positive correlation be- 
tween these two parameters (e.g. 3C273, 3C279, 3C454.4 and 
perhaps 0420—014). However, there are clear exceptions in 
either class, namely 0735+178 for the BL Lacs and 3C345 
for the quasars. 



5 DISCUSSION 

5.1 Overview: Shocks-in-Jets 

The most widely accepted model for the energy production 
and transport mechanisms in compact radio sources consid- 
ers the injection of material from the nucleus of the AGN 
into the base of the jet. As the plasma expands along the 
jet, flux-freezing implies that the component of the magnetic 
field perpendicular to the jet axis (B±) decreases with radial 
distance as r _1 , whereas the parallel component (By) should 
decrease more rapidly as r~ 2 (e.g. Laing 1981). Therefore, 
unless the field is being constantly re-isotropised, B± should 
become dominant. An enhancement of B± can also be pro- 
duced if an injection of relativistic particles into the jet leads 
to the emergence of a relativistic shock front. As the shock 
propagates down the jet, it energises the emitting particles, 
giving rise to flaring of the emission, and compresses the 
plasma, thus aligning magnetic field 'cells' with the direc- 
tion of the shock. The fact that we measure many more 
sources with B±_ than By could be evidence of these pro- 
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cesses taking place. On the other hand, B\\ could be en- 
hanced by shear at the edge of the jet caused by interaction 
with the surrounding medium (Laing 1981), or by the for- 
mation of highly oblique shocks in the jet (Cawthorne & 
Cobb 1990). However, the comparatively low levels of polar- 
isation measured at millimetre/submillimetre wavelengths, 
suggest that the magnetic field is still predominantly dis- 
ordered in the core regions of the jets. This could well be 
a consequence of the two orthogonal components interact- 
ing with each other due to the effects of these 'conflicting' 
physical processes. In most sources with intermediate lev- 
els of polarisation (3% < P < 9%) we have seen magnetic 
fields gradually changing direction. This allows us to specu- 
late that there may well be oblique shocks propagating along 
the jets, or that the apparent direction of the jets could be 
changing with time, or that the jets are curved on submil- 
liarcsecond scales. The best way to attempt to sort out the 
different possible causes of the apparent field rotation is a 
close examination of the multi-epoch data available on a 
number of sources. 

5.2 Plane Shocks 

If we start with the basic shock-in-jet scenario with no bend- 
ing, then for an initially disordered field, we predict that 
the formation of a shock results in an increase in the flux 
level, plus an increase in the polarisation with B± being en- 
hanced. Such behavour is exhibited (to some extent at least) 
by 0420-014, 0735+178, 1921-0293, 3C446 and 3C454.3, 
where the polarisation position angle tends to align itself 
with the jet when the flux is high, and deviates from align- 
ment as the flux decreases, but only 0735+178 and 3C454.3 
show a clear 90° swing in position angle and a positive cor- 
relation between flux and polarisation. On the other hand, 
most of the sources are not consistent with this simple pic- 
ture. In fact, 0235+164 shows the exact opposite behaviour, 
where the the magnetic field is perpendicular to the jet when 
the flux is low and becomes parallel as the flux increases. 
In this source, the flux and the polarisation are clearly 
anti-correlated. A few other sources (0528+134, 0736+017, 
1413+135 and 1749+096) show a similar tendency although 
the data available is scarce. A different pattern is observed 
in 3C279, BL Lacertae and 1335—127, whose magnetic fields 
are always perpendicular to the jet, despite the fact that all 
three have undergone very significant changes in flux and 
polarisation during the period covered by the observations. 
What is more puzzling, however, is that the relationship 
between P and F differs for each source: they are highly 
correlated in 3C279, anti-correlated in BL Lacertae and in 
1335—127 the correlation seems to change sign. 

If we modify the initial conditions so that the magentic 
field ordering prior to the formation of the shock is parallel 
to the jet (possibly due to shearing at the edge of the jet 
caused by interaction with gas in the host galaxy as sug- 
gested by Laing 1993), then we would predict an increase 
in the flux and a decrease in the polarisation as the shock- 
enhanced component of the magentic field B± cancelled out 
the polarised emission from the initially parallel field. The 
polarisation may then remain very low or could actually 
rise again, but this time with a position angle orthogonal 
to the initial direction. Thus, the overall behaviour could 
be complex with an initial anti-correlation between flux and 



polarisation followed by a period where they are correlated 
after a position angle swing of 90°. 

One of the best studied objects in our sample is 3C273 
and Stevens et al. (1997) combine some of our results with 
centimetre monitoring and VLBI to make a detailed compar- 
ison with the predictions of the shock model of Marscher & 
Gear (1985). The most significant observation about 3C273 
is that in late 1995 the flux and polarisation both increased 
as the source flared, however the initially parallel magnetic 
field remained parallel. Examination of the spectral index 
shows this cannot be due to an optical depth effect. This 
behaviour is not therefore consistent with the plane shock 
model predictions unless the source fortuitously undergoes 
a very large bend between the scales observed by centimetre 
VLBI and that of the millimetre/submillimetre emission. 

5.3 Oblique and Conical Shocks 

Lind & Blandford (1985) suggested that oblique and conical 
shocks could easily form in relativistic jets and showed that 
these geometries could have a dramatic effect on predicted 
number counts. These authors made little attempt to predict 
the direct effect of such shock structures on observations of 
individual objects. Cawthorne & Cobb (1990) however, in a 
prescient but largely- ignored paper, calculated the predicted 
effects of oblique and conical shock geometries on VLBI 
imaging and total polarised flux measurements of compact 
jet sources. They show that although in general such struc- 
tures will produce significant polarisation with position an- 
gle parellel to the jet (i.e. an enhancement of the perpen- 
dicular component of the magnetic field) they can also give 
rise to the opposite orientation for certain cone opening and 
viewing angles. A clear prediction of the model is that when 
this happens the polarisation will be low (< 10%). 

We believe that this type of shock geometry is entirely 
consistent with the behaviour we observe in many sources: 
generally low polarisation, often parallel inferred magnetic 
field and, furthermore, rapid apparent changes between or- 
thogonal magnetic field configurations which can be corre- 
lated or anti-correlated with increasing flux and/or polari- 
sation level. Examination of Figures 2 and 3 of Cawthorne 
& Cobb shows that small changes in either cone opening 
angle, or angle of the jet axis to our line-of-sight, can pro- 
duce a large swing in apparent magnetic field orientation. 
For example, in the case of BL Lacertae P and F are anti- 
correlated while the position angle changes only slightly 
about the jet's direction. OJ287 (another BL Lac object) 
exhibits the greatest variability in flux and polarisation in 
the whole sample. When the polarisation is high its position 
angle remains fairly unchanged at ~ 45° to the jet, despite 
large changes in flux. However, on three other occasions it 
flips from being parallel to perpendicular to the jet. 

5.4 Bending 

There are still sources in which it may be necessary to invoke 
bending of the jet. The best example is 3C345, whose jet is 
known to be very curved even in the highest-frequency VLBI 
maps. Our data covers the ocurrence of an outburst and its 
subsequent gradual decline to the pre-flare level. The degree 
of polarisation and position angle, however, show a lot more 
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variability (see Figure 3). This, together with the fact that 
the polarisation and the total intensity appear to be largely 
anti-correlated (see Figure 5) , is consistent with a jet being 
viewed at an initially small angle, which curves away from 
the line of sight as the shock evolves. This scenario is sup- 
ported by the prediction of Gopal-Krishna & Wiita (1992) 
that the percentage polarisation and total intensity are anti- 
correlated as the viewing angle increases within the range 
determined by the critical value 9 C = sirT 1 (1/r), where Y 
is the Lorentz factor. In fact, the jet of 3C345 has been es- 
timated to lie within 2° from our line of sight by Wardle et 
al. (1994) and as close as ~ 1° by Leppanen et al. (1995) 
(using the same analysis on a more recent VLBI map). The 
relatively low levels of polarisation measured by us on this 
source (< 4% on average) are also consistent with a small 
viewing angle. 

It is just possible that a conspiracy of bending on small 
scales might also reproduce the effects attributed to conical 
shocks in the previous section for OJ287, for example. This 
might then explain the earlier tentative result that P and F 
tend to be correlated in most quasars and anti-correlated in 
the BL Lacs as being due to the two groups tending to be 
viewed either side of the critical viewing angle 9 C . However, 
it is hard to explain the large degree of flux variability with 
small or no change in position angle, also commonly seen, 
in terms of bending. 

We suspect it is far more likely that whilst one or two 
sources do exhibit significant bending, in most cases the jet 
merely 'wiggles' a little about the axis, producing a modula- 
tion of the viewing angle of conical shock structures, whose 
opening angle also presumably evolves as it travels down an 
expanding jet. More theoretical work on the expected evo- 
lution of conical shock structures is probably required to 
make further more detailed comparisons of our results with 
the ideas of Cawthorne & Cobb. Stevens et al. (1997) have 
reached a very similar conclusion for the specific example of 
3C273. 



6 CONCLUSIONS 

(i) The magnetic field on sub-parsec scales in compact jet 
sources is less highly-ordered than on parsec scales, with 1.1 
mm polarisation levels less than 10% except in a very few 
sources, where we find it can be up to 15%. 

(ii) The magnetic field in the most highly-ordered sources 
is predominantly perpendicular to the jet axis, as measured 
on parsec scales. 

(iii) In the less-ordered sources, all possible apparent 
magnetic field orientations have been observed, as well as 
a high degree of variability of the measured polarisation po- 
sition angle in the multiply-sampled objects. These varia- 
tions are intrinsic and cannot be attributed to optical depth 
effects. 

(iv) In the best-sampled sources we find a good correla- 
tion or anti-correlation between total intensity and degree 
of polarisation. In less frequently sampled cases we also gen- 
erally find indication of such a correlation, and only in a few 
cases we find no correlation at all, although this may be due 
to insufficient sampling. 

(v) Our sample includes similar numbers of BL Lacs and 
high-polarisation quasars, as well as a few low-polarisation 



quasars and radio galaxies. We find, statistically, no differ- 
ence between the degree of ordering or magnetic field geom- 
etry of the BL Lacs and high-polarisation quasars. 

(vi) Although the behaviour of a few sources (again 
mainly the highly-ordered cases) is very consistent with the 
predictions of a plane shock-in-jet model such as that of 
Marscher & Gear (1985), in most cases the behaviour is in- 
consistent with this picture, unless very large bends occur 
on scales smaller than those so far sampled by VLBI. 

(vii) The behaviour of these sources can be explained by 
the conical shock model of Cawthorne & Cobb (1990). How- 
ever it may still be necessary to invoke some degree of bend- 
ing in certain cases. 

We have shown the power of a new technique, that 
of millimetre/submillimetre polarimetry, in deducing new 
constraints on models of compact, flat-spectrum, radio-jet 
sources, which have been well-studied using other techniques 
for many years. One of the biggest limitations on this ex- 
periment has been the uncertainty in the jet direction on 
the very small (milliparsec) scales sampled at these wave- 
lengths. Polarisation-sensitive VLBI measurements at 1 mm 
would be the ultimate tool for this experiemnt, however that 
is likely to be some years away. A combination of the high- 
est possible frequency VLBI and millimetre-submillimetre 
photopolarimetric monitoring with further theoretical de- 
velopment is likely to be the way forward in the near future. 
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APPENDIX A: SOURCE SUMMARY 

0133+476 (OC'457) 

A single observation is available on this source giving a flux 
density of 1.5 Jy and a percentage polarisation of 10.5%. 
This is a 6a measurement obtained in only 25 minutes inte- 
gration time on a relatively faint source. For a 9y lb i = 165° 
obtained from a low resolution 5 GHz VLBI map, the in- 
ferred magnetic field orientation is 18° away from the per- 
pendicular to the jet. 

0235+164 (OD160) 

There are data on this source for four epochs obtained in 
just over a year. The flux density changes from about 1 
to more than 3 Jy while the polarisation drops from 6 to 
less than 3%, indicating a possible anti-correlation between 
these two parameters (r = —0.6 at 60% significance level). 
The mean flux and polarisation over the four epochs arc 
2.2 Jy and 4.9% respectively. A three-point spectral index 
has been obtained on three epochs giving a shallower than 
average spectrum (relative to the other sources in the sam- 
ple). The position angle flips by about 90° at the epoch of 
highest flux. This is a very compact source with no clear 
milliarcsecond structure, but for a tentative jet orientation 
Ovlbi = —150° obtained at 22 GHz, the magnetic field 
would appear to be perpendicular to the jet during the first 
two epochs, becoming aligned with it after the flare. 

0316+413 (3C84) 

This source is a rather peculiar radio loud seyfert. Although 
it is fairly bright (of the order of 2 Jy at 1.1 mm) it seems 
to have very low polarisation in comparison with most of 
the sources in the present sample (~ 1.5%). This makes it 
very difficult to achieve polarimetric observations with good 
signal to noise. Only two detections have been obtained ten 
months apart to just over the 2a level, which do not show 
any significant changes in the emission during that period. 
High-frequency VLBI maps are available which indicate an 
inner jet position angle of —150°; the inferred magnetic field 
orientation is neither perpendicular nor parallel to the jet. 

0415+379 (3C111) 

Although a ~ 10<r detection was obtained on this source at 
all three photometric wavelengths and the source is fairly 
bright (>2 Jy at 1.1 mm) only a la upper-limmit was mea- 
sured in polarimetry at the 0.7% level. 

0420-014 (OF- 035) 

This source was observed four times. The mean polarisa- 
tion and flux density are 5% and 2.8 Jy. One year separates 
epochs 2 and 4, during which the 1.1 mm flux drops steadily 
to almost a quarter. This does not reflect greatly on the po- 
larisation level but it rotates the position angle by more than 
50°. The core of this source is largely unresolved but very 
recent 43 GHz VLBI images suggest a jet position angle of 
— 150°. This indicates that the magnetic field tends to align 
itself with the jet as the flare loses strength. 

0430+052 (3C120) 

The flux density measured on this radio galaxy at 1.1 mm 
is 0.75 Jy and only a la upper-limmit of 2.2% was obtained 
in polarimetry. 
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0528+134 (OG147) 

There are two epochs of observations on this source, taken 
four months apart, during which the polarisation drops by 
40 per cent as the flux almost doubles (increasing from 1.9 
to 3.5 Jy) and the position angle suffers a rotation of 80°. 
VLB1 observations at 22 GHz suggest a jet position angle 
close to 70° which indicates that the magnetic field turns 
towards alignement with the jet as the flux increases. 

0735+178 (01158) 

There are four epochs of data for this source. The mean po- 
larisation and flux density are 6.3% and 1.2 Jy. The milliarc- 
second jet orientation is reasonably well determined from 
VLBI data (Ovlbi = 73°) and the magnetic field is seen to 
change from being perpendicular to being parallel to the jet 
as the flux decreases from ~ 1.5 to 0.8 Jy. The polarisation 
also drops during the same period. 

0736+017 (OI061) 

This is the faintest source on which we have obtained po- 
larimetric detections: with a flux density of < 0.6 Jy the 
polarisation measured was 5.4% (with signal to noise > 2a). 
An increase in flux to 0.76 Jy four months later is accompa- 
nied with an increase in percentage polarisation to 14.2% (a 
7<r detection). Between these two epochs the position angle 
is rotated by 43° . The jet position angle observed on a re- 
cent 22 GHz map (Ovlbi = —80°) suggests that as the flux 
and polarisation increase, the magnetic field tends towards 
alignement with the jet. 

0829+046 

This source has been measured twice within four months, 
during which the flux density increased from 0.93 to 1.44 Jy 
and the percentage polarisation decreased from 4.9 to 2.9%, 
i.e. they seem to be anti-correlated. The position angle re- 
mained largely the same but the orientation of the magnetic 
field cannot be infered since no VLBI observations of this 
source have been found in the literature. 

0851+202 (OJ287) 

This source was observed on eight different occasions and 
in terms of its polarisation properties it is the most variable 
source in the sample. The polarisation ranges widely be- 
tween 5.1 and > 15% averaging to 9.8% over all epochs. 
However, its measured flux density only varied between 
1 and 2.4 Jy. The milliarcsecond jet position angle has 
been determined from very-high-frequency VLBI maps to 
be Ovlbi = —120°. Over most epochs, the magnetic field 
geometry inferred from the measurements of the polarisa- 
tion position angle is neither parallel nor perpendicular to 
the jet. During the last two epochs, however, the polarisa- 
tion drops from 13% to 5% while the flux increases from 
1.3 to 2.4 Jy. The position angle is clearly seen to rotate 
by roughly 90°, suggesting that the magnetic field changes 
from being perpendicular to the jet (at the high polarisation 
epoch) to being parallel to it. Over all epochs no correlation 
is found between the changes in flux and degree of polarisa- 
tion, but if the data are examined between epochs 1—5 and 
5—8 separately, there is a hint of a positive correlation first 
r = 0.7 and an anti-correlation later r — —0.8 both with the 
same 80% confidence level. 



0923+392 (4C39.25) 

There are four polarisation and five photometric epochs on 
this source. The variation in polarisation (and flux) is not as 
dramatic as seen in other sources, ranging between 3.7 and 
5.6% with measured position angles spread over a ~ 45° 
range. However, a flare was recorded between the last two 
epochs (taken one month apart) when the flux rose from 
a minimum of 1.1 Jy to a maximum of 2.5 Jy. The core's 
position is uncertain since it remains unresolved with VLBI. 
Relative to a position of —88° for the innermost resolved 
component, the magnetic field is neither perpendicular nor 
parallel to the jet, except just before the flare where it is 
found within 15° of the aligned orientation. 

1055+018 (4C01.28) 

Two observations are available on this source, taken almost 
three years apart, which are very similar: the polarisation is 
very high (9% on average) and the mean flux is ~ 1.5 Jy. 
The magnetic field appears to be exactly perpendicular to 
the jet, however the jet position angle (305°) was obtained 
from low- frequency maps (1.7 GHz) and it may not represent 
the true orientation of the milliarcsecond jet. 

1226+023 (3C273) 

This source was observed at ten different epochs. The flux 
density at 1.1 mm varied in the range 11 - 21.2 Jy. This 
unusually high brightness allows polarisation levels as low 
as 0.7% to be measured with enough signal to noise, so far 
not achieved on any other source. In the last two epochs (in 
late 1995) this source was caught flaring and the polarisa- 
tion rose to 6.8%. Except at the epoch of lowest polarisation 
(and flux) the position angle lies within a small range of 30° . 
From numerous high-frequency VLBI observations, we esti- 
mate the orientation of the milliarsecond jet to be ~ —110°, 
however at 100 GHz the inner jet exhibits a somewhat wiggly 
morphology. The inferred magnetic field geometry tends to- 
wards alignment with the jet, except at the low polarisation 
epoch when it is closer to the perpendicular configuration 
(although with very low level of ordering). There is also a 
significant correlation between flux and degree of polarisa- 
tion with r — 0.8 at the 99.4% significance level. 

1253-055 (3C279) 

This is the most frequently observed source in the sample 
with 11 polarisation epochs. It is very bright (with flux den- 
sity in the range 6 - 13.5 Jy) and very highly polarised (with 
an average of 11.3%). High-frequency VLBI maps reveal an 
inner jet position angle of —135° within O.lmas. The position 
angle varies only by a few degrees between epochs, indicating 
a highly-ordered magnetic field unambiguously perpendicu- 
lar to the jet. The flux and polarisation are highly correlated 
with r = 0.76 at the 99.4% significance level. 

1308+326 

A 4<r observation was obtained on this source giving a per- 
centage polarisation of ~ 10% and a flux density of 1 Jy. Us- 
ing a VLBI jet position angle of —46° (obtained at 5 GHz) 
the inferred magnetic field orientation is neither parallel nor 
perpendicular. However, the jet curves significantly within 
the inner 4mas and it is likely to bend further within the 
unresolved region, casting some doubt on the validity of the 
Ovlbi used. 
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1335-127 (1334-127) 

There are five observations on this source with a flux density 
in the range 2.5 - 4.9 Jy, which overall appears to be anti- 
correlated with the polarisation that ranges between 2.5 and 
7.9%. The position angle is very steady remaining within a 
23° range over all epochs. A very tentative jet orientation of 
135° has been obtained from a 5 GHz VLBI map in which 
the core is not clearly resolved. This would indicate a mag- 
netic field geometry perpendicular to jet. 

1413+135 

This source has been observed three times, the last two 
epochs taken within a month. The flux density ranges be- 
tween f .1 and 1.6 Jy and appears to be anti-correlated with 
the polarisation that varies in the range 3.7 - 9.1%. The 
only VLBI map available is of low resolution but shows a 
double-sided jet, the most compact side of which is oriented 
along —100° at ~ 8mas. This is probably a poor estimate 
but the indication is that the magnetic field geometry lies 
roughly parallel to the jet when the polarisation is high. 

1514-241 (AP Lib) 

This is a low declination source not frequently observed. The 
only observation obtained gives a flux density of 1.1 Jy and 
polarisation of 5.5% to better than 3a. The position angle 
of the jet has only been estimated from a VLA map with 
insufficient resolution close to the core to try to establish the 
orientation of the magnetic field. This measurement gives a 
very flat spectral index, in broad agreement with the shal- 
lower than average spectral index values obtained on this 
source by Gear et al. (1994) over three different epochs. 

1633+382 (4C38.41) 

The single observation available on this source was obtained 
with a flux density <0.9 Jy, however the polarisation was 
high at ~ 8% with a signal to noise of 3a. The polarisation 
position angle indicates a possibly perpendicular magnetic 
field relative to a jet orientation of —45° determined from a 
low-frequency (5 GHz) VLBI map. 

1641+399 (3C345) 

This source has been observed eight times. The degree of 
polarisation changes readily between the extremes of 0.7 and 
11.2%, averaging to 4.7% over all epochs. The flux variability 
is very different, since a large flare at the beginning of the 
observations (when the flux increases from 1.6 to 5.8 Jy) 
is followed by a slow decay to its pre-flare level over the 
subsequent three years. The jet is known to be very twisted 
and it is not possible to determine its orientation close to the 
core with any certainty. High-frequency VLBI images show 
that emerging components follow paths in the range —135 
—60°. The polarisation position angle changes by tens of 
degrees between epochs. Using a recent value for Ovlbi = 
—55° (obtained at very high frequency) the magnetic field 
is found within 30° of alignment with the jet on six epochs, 
but clearly perpendicular to it on one occasion. There is a 
90° difference in position angle between the first and last 
polarisation epochs, when the polarisation was close to 5% 
in both cases but the flux measured was 5.8 and 1.9 Jy 
respectively. 



1741-038 (OT-68) 

There is a single observation on this source with just under 
3a signal to noise. The polarisation is 4.6% with a position 
angle of 126° and the flux density is 1.16 Jy. The inferred 
magnetic field orientation is parallel to the jet with 8vlbi — 
-155° (recently observed with VLBI at 22 GHz). 

1749+096 (OT081) 

This source was observed on three occasions over 18 months. 
The polarisation rises from 3.7 to 8.5% while the flux drops 
from 4.6 to 1.9 Jy, i.e. they are anti-correlated. The posi- 
tion angle varies within a 50° range, but the compactness 
of the core makes the jet position angle 65° used somewhat 
uncertain and the magnetic field geometry difficult to infer. 

1823+568 (4C56.27) 

There is only one epoch on this source which is highly po- 
larised at 12.8% and with a flux density of 1.1 Jy. The mag- 
netic field seems to lie within 25° from the perpendicular to 
a jet orientation of —166° determined from 5 GHz maps. 

1921-293 (OV-236) 

This low declination source was observed on four different 
epochs. It is very bright, with a flux density ranging from 

4.7 to 7 Jy, and highly polarised with a mean polarisation 
level of 9.1%. The position angle changes by tens of degrees 
between epochs but the magnetic field remains largely per- 
pendicular to the jet (only one VLBI map is available but it 
has been obtained very recently at 43 GHz). The flux and 
degree of polarisation are uncorrelated. 

2200+420 (BL Lacertae) 

BL Lacertae has been observed six times. It is not brighter 
than the other BL Lacs with a mean flux of 1.8 Jy (varying 
between 1 and 2.6 Jy) but it is the most highly polarised 
with a mean polarisation of 8.1%. The changes in flux and 
percentage polarisation are anti-correlated with r = —0.77 
at 93% confidence level. The polarisation position angle is 
the same on four of the six epochs, thus remaining largely 
within a level consistent with a magnetic field geometry per- 
pendicular to the jet. The jet orientation is well determined 
from several high-frequency VLBI maps. 

2223-052 (3C446) 

There are two epochs on this source, taken more than two 
years apart, between which the flux density and the polar- 
isation both drop significantly from 4.5 to 1.1 Jy and from 

3.8 to 2.2% respectively. The position angle of the jet is 218° 
(from a high resolution 100 GHz VLBI map) and suggests 
a perpendicular magnetic field when the source flux is high, 
deviating from this position by 50° when the emission is low. 

2251+158 (3C454.3) 

This source has been measured six times. The flux and po- 
larisation lie in the ranges 3.8 - 6.7 Jy and 2.9 - 8.3% re- 
spectively, and are very strongly correlated (r = 0.94 with a 
99.5%confidence level). The position angle changes readily 
within a > 90° range. The epoch of lowest flux and minimum 
polarisation coincide with a magnetic field aligned with the 
jet, which has a position angle of —65° (from 43 GHz VLBI 
maps). As the flux and polarisation increase the magnetic 
field turns towards the perpendicular configuration. 
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